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High Multiplicity pp and pA Collisions are a place where the macroscopic description (thermo- 
dynamics and hydrodynamics) meets with the microscopic one (pomerons and QCD strings). In 
the first part of the paper we study what happens with the hydrodynamical predictions as the 
system size gets smaller and smaller. For simplicity, we don't do it numerically, but analytically 
using Gubser' s flow. We found that the radial flow is expected to increase, while the elliptic flow 
decreases, and high harmonics become perhaps too small to be observed. In the second part of the 
paper we approach the problem from the opposite side, using a string-based Pomeron model. We 
found that as the intrinsic temperature of the string grows, it approaches the Hagedorn regime and 
produces a high entropy string ball, amusingly dual to a certain black hole. Furthermore, when the 
string temperature narrows on the Hagedorn temperature or T/Th — 1 = 0(1/N C ), the stringy ball 
becomes a sQGP ball with non-negligible pressure and hydrodynamical flow. 



I. INTRODUCTION 
A. The outline 

High energy heavy ion collisions are theo- 
retically treated very differently from pp and 
pA ones. While the former are very well de- 
scribed using macroscopic theories - thermo- 
dynamics and relativistic hydrodynamics - the 
latter are subject to what we would like to call 
the "pomeron physics", described with a help 
of microscopic dynamics in terms of (ladders 
of) perturbative gluons, classical random gauge 
fields, or strings. The temperature and entropy 
play a central role in the former case, and are 
not even mentioned or defined in the latter case. 

The subject of this paper is the situation 
when these two distinct worlds (perhaps) meet. 
In short, the main statement of this paper 
is that specially triggered fluctuations of the 
pp and pA collisions of particular magnitude 
should be able to reach conditions in which the 
macroscopic description will also become pos- 
sible. While triggered by experimental hints at 
LHC to be discussed below, this phenomenon 
has not yet been a subject of a systematic study 
experimentally or theoretically, and is of course 
far from being understood. Needless to say, the 
problem touches on a number of key issues, we 
now enumerate and then address. From the 



macroscopic side, they are: 

• What is the smallest system size which 
still undergoes a hydrodynamical explo- 
sion? 

• How do all hydrodynamical observables 
scale with the system size R and viscosity- 
to-entropy ratio r]/s, for such systems? In 
particular, how large are the viscous cor- 
rections for radial and elliptic flows? 

• How do amplitudes of higher angular har- 
monics v n scale with n,R and r\jsl In 
which p t region do we expect hydrody- 
namics to work, and for each v n l 

• Do high multiplicity pp and pA collisions 
in which the (double) "ridge" has been re- 
cently observed at LHC [7H1] fit into such 
a hydrodynamical scaling? 

From the microscopic side we will use a 
stringy description of the pomeron at strong 
coupling in the context of holographic QCD 
with a confining wall. Specifically, we will use 
the recent model developed by Staffers and Za- 
hed [TJ |2] where both the issue of a coherent 
scattering amplitude and entropy production 
can be adressed simultaneously. The issues we 
will address in this model are: 
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• What is the "intrinsic" string temper- 
ature, for the strings produced in the 
minimal-bias (unconstrained) pp and pA 
collisions? How does the multiplicity de- 
pend on the collision energy for the usual 
pomeron? 

• What difference in the string shape and 
dynamics should be observed when the 
string becomes "near-critical" as its in- 
trinsic temperature T approaches the 
Hagedorn temperature Tjj? 

• What are the energy, pressure, entropy 
and size of such a "near-critical" string? 

• Is the thermodynamics of the near-critical 
string "holographic" ? (Is it dual to a cer- 
tain black hole?) 

• If one uses the Kubo formula to evaluate 
the "string viscosity" , what values are ob- 
tained ? 

• For extremely strong fluctuations, with 
T - T H so small that the "string ball" 
thermodynamics changes from the Hage- 
dron regime to the partonic (sQGP) one, 
what behavioral change is expected? 

B. Phenomenology 

The major objective of the heavy ion col- 
lision program is to create and study proper- 
ties of a new form of matter, the Quark-Gluon 
Plasma. Among the many proposed signatures 
priginally in |3], a prominent place is taken 
by a collective "explosion", described by rel- 
ativistic hydrodynamics. A major discovery 
at RHIC has been a robust radial and elliptic 
flow, in excellent agreement with hydrodynam- 
ics. This point is now also supported by LHC. 
The flow is characterized by higher angular mo- 
ments v n (pt) in the range of p t < 3GeV, also 
in good agreement with hydrodynamical pre- 
dictions. 

One may naturally ask if large size colliding 
nuclei (AuAu at RHIC and PbPb at LHC) are 
really necessary for this description to hold, or 
if one can achieve this for smaller systems as 
well, such as those in pA or even pp collisions. 
In fact radial flow effects were searched for more 



than 30 years ago, in the minimum-bias pp col- 
lisions at CERN ISR [I], with negative results. 
Some weak indications for the radial flow have 
been found in specially triggered pp collisions 
by the FERMILAB MINIMAX experiment [5], 
but the data remained inconclusive. 

The "ridge" itself is an angular correlation 
in the azimuthal angle between two particles at 
Acf> < I which extends to large rapidity range 
|Ay| > 4. Well known in AA collisions, it has 
been discovered by the CMS collaboration [5] in 
the hrst LHC run. In order to see it, one has to 
apply a very restrictive rare trigger. The CMS 
pp study had selected out of ~ 10 9 events only 
~ I0 3 , populating the highest multiplicity bin 
with a convincing ridge effect. More recently 
the same phenomenon was seen in pPb colli- 
sions as well, now by the CMS [7J, ALICE [S] 
and ATLAS [9 collaborations. Larger number 
of "participant nucleons" and higher average 
multiplicity substantially weaken the require- 
ments for the trigger: the "ridge" is seen at the 
trigger level of few percents higher multiplicity 
events. It is shown in those works that in pp 
and pA collisions, the same threshold in terms 
of multiplicity is needed to start showing the 
"ridge" . 

Angular correlations naturally appear as a 
consequence of a hydrodynamical explosion of 
a non-azimuthally symmetric objects. The spa- 
tial shape is then translated to momentum 
space and is observed. For example, in the com- 
ments on the CMS discovery written by one of 
us [10] it was illustrated by a string placed out- 
side of an (axially symmetric) "stick of explo- 
sive. While the basic wind blowing is isotropic 
in <f), an extra string may move in a preferred 
direction. In central AA collisions it is similar 
to that. A symmetric explosion has perturba- 
tions in the form of localized "hot spots" . But 
in general, any sufficiently deformed initial col- 
lisions for the fireball would be sufficient to cre- 
ate ridge-like correlations. 

Furthermore, the subtraction of the so called 
"back-to back recoil" (a peak at <j> ~ ir) which is 
either evaluated from some mini jet model (HI- 
JING) or studied experimentally in the smaller 
multiplicity bins, reveals that a ridge is "dou- 
bled. The remaining correlation function is 
nearly symmetric with and tf> w it show- 

ing nearly the same correlation. In other words, 
there is a clearly dominance of the second angu- 
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lar harmonics [8, 9 . (This is unlike the central 
AA data, in which the strongest harmonics is 
the third n = 3.) 

A number of authors described this phe- 
nomenon hydrodynamically. For the pA case, 
one may mention e.g. which starts from 

Glauber-inspired initial conditions similarly to 
what is done in the AA case. Indeed, one pro- 
ton propagating through the diameter of the 
Pb nucleus "wounds" 10-20 nucleons. Those 
are placed in certain random pattern in the 
transverse plane, providing some estimates of 
the higher angular harmonics e n from which hy- 
drodynamical perturbations may start. The re- 
sults are qualitatively consistent with the LHC 
data. 

The objective of this paper is to extend such 
studies, using instead of a complicated "real- 
istic models" with huge number of details and 
heavy numerics (the "event-by-event" hydrody- 
namics) an analytic approach. As we will see, 
this will allow us to focus on generic depen- 
dences of the predictions on the parameters of 
the problem. 



C. Hydrodynamics 

Before going to the details, let us express the 
main ideas and tools we have in mind in general 
terms. Since we will be discussing the valid- 
ity of hydrodynamics, it is important to specify 
from the onset which hydrodynamics we will 
use. 

The simplest are the equation of ideal hydro- 
dynamics 



d^ v = 



(1) 



containing the stress tensor in local approxima- 
tion (no derivatives). Thus the only deriva- 
tives appear linearly and therefore simultane- 
ous rescaling of the size and the evolution time 
x^ — > Xx^ does not change the equation. So, 
ideal hydrodynamics will produce the same so- 
lution for fireball of any size, provided other 
parameters are unchanged. And indeed, such 
scaling was, for example, seen experimentally, 
e.g. for AuAu and CuCu collisions at RHIC. 
However, this works only if in both cases all 
viscous corrections are negligible. 

At the very least, this requires as one applica- 
bility condition that the system's size R is much 



larger than some microscopic scales such as the 
inverse temperature T _1 . The corresponding 
ratio will be one small parameter which will 
subsequently change 



0(1) 



f jPJ£pP,high mutiplicity 

l l 

> „„ t — r > 



^central pA J^J^central AA 



(2) 

0(1/10) 



we 



Another important small parameter which wc 
seem to have for strongly coupled Quark-Gluon 
Plasma (sQGP) is the viscosity-to-entropy- 
density ratio 



0(1/10) < 1 



(3) 



Roughly speaking, it tells us that viscous ef- 
fects - or the mean free path - is additionally 
suppressed compared to the micro scale 1/T. 
The smallness of this second parameter sug- 
gests that one can hope to apply hydrodynam- 
ics even for systems as small as pp collisions, for 
which the former parameter is no longer small. 

Hydrodynamical effects can be further sep- 
arated into an overall explosion - called the 
radial flow, and higher angular harmonics of 
flow. The former is expected to be enhanced 
by larger gradients, while the latters get more 
suppressed by viscosity effects. Finally, one can 
discuss the applicability of hydrodynamics as a 
function of the particle transverse momentum 
v„(p t ). 

We will show below that while in central 
AA we see about 6 angular harmonics, in pA 
only the second one remains relatively weakly 
damped by viscosity, while in the pp case even 
the second one is rather strongly damped. We 
still expect the radial flow (zeroth harmonics) 
to be enhanced, as the bulk viscosity is not sup- 
posed to be large enough in this case. 

Furthermore, we will argue that the high 
multiplicity events in pp should be attributed to 
a QCD string going near- critical, with T — > Th- 
If so, the transverse size R p f is of a sub-nucleon 
scale and is related to the significantly smaller 
QCD string scale 



'a' 



(2GeV)~ 



0.1 fm 



(4) 
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D. Microscopic models 

The usual (min-bias) pp and pA collisions 
are described by a number of models which use 
pomeron and reggeon exchange. Those terms 
originated from certain parameterization of the 
elastic and diffractive scattering amplitudes at 
high energy kinematics \fs 3> \f—i. The micro- 
scopic derivation in weak coupling is done via 
a (rapidity ordered) BFKL gluon ladders [12]. 

At strong coupling it is dominated by some 
form of a string (surface) exchange. In the 
AdS/CFT holographic context the discussion 
of elastic and inelastic amplitudes generated 
by such surface exchanges were originally ad- 
dressed using bosonic variational surfaces [13] 
[14] . It has also been realized that the holo- 
graphic pomeron can be associated with a spin- 
2 graviton exchange [J5] , in pure AdS with su- 
persymmetry and conformal symmetry, which 
has triggered a rather extensive literature. In 
particular, a black-disk model [15] . 

Following on the semi-classical surface ex- 
changes in [T31 [H] , a quantum bosonic string 
(surface) exchange model was suggested re- 
cently in Witten' s confining background [T7] . 
We will use its variant as discussed in [J [5] in 
AdSs with a confining wall where the impor- 
tant number of transverse directions is physi- 
cally identified. We will refer to this as the SZ 
model. Let us preview its physics in some gen- 
eral terms and detail it later. 

One important idea is that the string which 
dominates semiclassically the scattering ampli- 
tude has some excitations, in the one-loop ap- 
proximation. The quantum action helps to 
explain why the pomeron is supercritical and 
provides a specific assignment for its intercept 
ap(t = 0) — 1. These excitations of the string 
modes can be described by a one parameter, an 
intrinsic (Unruh) temperature T. This temper- 
ature is set by the kinematics of the minimal 
bias collisions, but it can also fluctuate describ- 
ing less usual events. When this temperature 
is small as compared to the so called Hagedorn 
temperature T <C Th, the string fluctuations 
are small and can be treated as a set of oscilla- 
tors, as we schematically show in FigJI^. This 
is the case for the traditional applications men- 
tioned above, to the minimum-bias pp collisions 
at high energies up to LHC. 



a b 

FIG. 1: String exchange between two sources 
(crosses) separated by the impact parameter b: the 
cold string case j3 < /3h (a); the near-critical string 
case /3 —> Ph (b) . 



With increasing collision energy the effective 
temperature grows and for some super-high col- 
lision energies (not reached at colliders) it may 
approach the Hagedorn temperature T — > T^. 
At current energies (LHC) it can also happen, 
as fluctuations. We will argue that in this new 
regime the string will develop large excitations 
in the form of a "string ball" depicted in FigJTJj. 

This near-Hagedorn regime is well known 
in finite-temperature QCD as the Polyakov- 
Susskind phenomenon and the coresponding 
entropy behavior is schematically shown in 
Figj2] The entropy and energy asymptote in- 
finity at T 4 Th due to the proliferation of 
string states. The free energy and pressure re- 
main finite. The rapid rise will be referred to 
as the "near-critical" or Hagedorn regime, also 
known as an interval with a "soft" equation of 
state 

T « T H , p«e (5) 

For a general review discussing effective strings 
and their role in the gluodynamics phase tran- 
sition at finite and large N c see e.g. [TH]. 

Since in the prompt collision the pomerons 
override the reggeons at large y/s, the colli- 
sion is dominated by gluons instead of quarks. 
Therefore the pertinent Hagedorn regime is 
that of gluodynamics rather than QCD. We re- 
call that the critical temperature in gluodynam- 
ics is T c yM w 270 MeV, which is significantly 
higher than the critical temperature in QCD 
which is T® CD « 165 MeV. The former is im- 
portant at the formation stage while the lat- 
ter is important at freeze-out. In the effective 
string language, the transition in gluodynamics 
involves closed strings while that in QCD in- 
volves open strings. If we recall that the ratio of 
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the closed to open string tensions is gaI^f ~ 2 
and that the Hagedorn temperature is propor- 
tional to the square root of the string tension, 
we expect T? M /TQ CD ps y/2, which is about 
the experimentally measured ratio. 

We suggest that the behavior of a "string 
ball" , created on a pomeron string as a fluctu- 
ation, near the Hagedorn temperature or T ps 
Th, should be very similar to that of strings 
in thermal equilibrium and close to deconfine- 
ment. The latter is referred to as the mixed 
phase of pure gauge theories with N c > 2. In- 
deed, the gluodynamics transition is dominated 
by a dilute gas of close strings each of which 
carry a considerable stringy entropy. At the 
Hagedorn point, this stringy entropy is best 
packaged in a large and space filling string much 
like our string ball. 

We further suggest that this Hagedorn phe- 
nomenon explains the significant increase of the 
produced entropy and thus multiplicity of sec- 
ondaries. We believe its possible production 
in high-multiplicity pp collision is a new phe- 
nomenon not discussed previously. Another 
new element of our discussion (which is based 
on some recent ideas and technical progress 
in string theory) is the strong similarities we 
demonstrate between this "string ball" and the 
black hole, in terms of an effective temperature- 
entropy relations and even an effective viscosity 
we will evaluate. 

Finally, we will argue that when the differ- 
ence between T and Tjj becomes very small and 
the energy (entropy) densities become as large 
as 



s 

J-3 



T C T H 



FIG. 2: Schematic temperature dependence of the 
entropy density. The dashed line represents equi- 
librium gluodynamics with a first order transition 
at T = T c . The solid line between points A and B 
represents the expected behavior of a single string 
approaching its Hagedorn temperature Th- 



perhaps the first observed experimental mani- 
festation of this phenomenon. 



II. HYDRODYNAMICS AT ITS EDGE: 
STUDIES IN THE NAVIER-STOKES 
APPROXIMATION 



Ideal radial flow using Gubser's 
solution 



1 H 



T 3 

1 H 



Nz 



(6) 



a second qualitative change happens: the 
"string ball" becomes explosive, as the pressure 
is no longer small. We argue by analogy with 
the thermodynamics of gauge theories: at some 
entropy density the mixed phase ends and the 
deconfined phase takes over. While we know it 
is so in equilibrium, we suggest the same should 
happen out of equilibrium, in a small fireball 
originating from a string ball in a Pomeron. Its 
direct observable consequence should be the ap- 
pearance of a significant pressure p ps e/3 and 
subsequent hydro explosion. The appearance 
of the (double) ridge in pp and pA collisions is 



Many groups have studied heavy ion colli- 
sions by solving hydrodynamical equations nu- 
merically. However, it is not really necessary 
for the current purposes since there exist a rel- 
atively simple analytic solution found by Gub- 
ser [T5], see also [5U]. In this section we will 
apply this solution and compare the AA hydro- 
dynamics with pA and pp ones. 

Gubser flow is a solution which keeps the 
boost-invariance and the axial symmetry in the 
transverse plane. It is obtained via symmetry 
under special conformal transformation, and 
therefore, the matter is required to be confor- 
mal, with the EOS 



e = 3p = T /* 



(7) 



G 



where the parameter /» = 11 is fitted to repro- 
duce the lattice data on QGP thermodynamics. 

The coordinates used are either the usual 
proper time -spatial rapidity - transverse radius 
- azimuthal angle set (f , 77, r , <p) with the metric 

ds 2 = -d? 2 +? 2 dri 2 + df 2 + f 2 dc/) 2 , (8) 

or the comoving coordinates we will introduce a 
bit later. The solution will depend only on the 
dimcnsionful variables f , f. The bar disappears 
because the solution has only one parameter q 
with dimension of the inverse length, so one can 
use rescaling to dimensionless variables 



t = qr, r = qr 



(9) 



Gubser's solution of ideal relativistic hydro- 
dynamics, for the transverse velocity and the 
energy density reads 



vj_(t,r) = 



2tr 



l + t 2 



e 2 8 / 3 



(10) 



9 4 t 4 /3[l + 2(<2 +r 2) + ( <2 _ r 2)2]4/3 

where eo is some normalization parameter. 

The first task we perform is mapping the AA, 
pA and pp collisions to these coordinates, which 
is done via selection of the scale factors chosen 



as 



Qaa = 4-3, g-j = 1, g" 1 - 0.5 (fm) (11) 

and the energy density parameter, which is re- 
lated to the entropy-per-rapidity density of the 
solution 



(a 



= f- 1 / 3 ( 3 dS \ 
\ 167r dr\ ) 



4/3 



(12) 



As in the Bjorken scaling solution, the rapidity 
interval is infinite and nothing depends on 77, 
so the total entropy is infinite. Rapidity inde- 
pendence is well satisfied in the LHC data we 
discuss, in the observed rapidity interval. The 
entropy density is directly related to observed 
density of charged secondaries via a relation 



dS 
drj 



7.5 



dN ct 
drj 



(13) 



defined at freezeout. 

We use for central LHC AA=PbPb collisions 



dNf h A /di] = 1450 



(14) 



(Note that we ignore the difference between 
rapidity and pseudo rapidity.) The pp and 
pA data are split into several multiplicity bins: 
for definitcness, we will refer to one of them 
in the CMS set, with the (corrected average) 
multiplicity N c h = 114 inside \r)\ < 2.4 and 
Pt > 0.4 GeV acceptance. We thus take 

dNtf/dr, = dN%/d V = 1.6 (15) 

where the factor 1.6 approximately corrects for 
the unobserved p t < 0.4 GeV region. 

Now the energy parameters are also fixed, 
and one can map the region in which hydro- 
dynamics is supposed to work for all three col- 
lision regimes. Such a map is shown on the t, r 
plot in Figj3] Hydrodynamics is assumed to be 
valid between the (horizontal) initial time lines 
and the contours of fixed freeze out tempera- 
ture TV. Note, that while the absolute sizes 
and multiplicity in central AA are quite dif- 
ferent from pA, in the dimensionless variables 
those are closer to each other than in pp. The 
reason is an assumed factor of 2 jump in the 
radius. 

Thus the cases get more and more "explo- 
sive" , as the gradients progressively grow. The 
spectra should be calculated by the standard 
Cooper-Fry formula 



dN 



d-qdp 



p^dY;^ exp 



(16) 



in which is the freeze out surface, on which 
the collective velocity u^(t, r) should be taken, 
for details see 124 . The transverse collective 



velocity on the freeze out curve is read off ( 10 1 
which is too simple to plot. We just mention 
the maximal values reached at the "knee" of 
the freeze-out curves. They are 



v^ ax [AA lP A,pp] 



[0.69,0.83,0.95] (17) 



respectively. For large pt of the secondary a 
region around this "knee" dominates the in- 
tegral. The slopes of the spectra of the rela- 
tivistic pion are simply blue-shifted. Spectra 
of heavier particles - protons and antiprotons 
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becomes an ordinary differential equation 



,dT 



FIG. 3: (color online) The three horizontal lines 
correspond to the initial time: from bottom up AA 
(blue solid), pA (dash black) and pp (red dash-dot). 
The corresponding three curves with the same color 
are the lines at which the temperature reaches the 
same freeze-out value, set to be Tf — 150 MeV. 
The two thin solid lines correspond to the values 
of the variable p = — 2.2 (lower) and —0.2 (upper). 
Those values are used as initial and final values in 
the evolution of higher harmonics. 



especially - are modified differently, leading to 
a strong enhancement of the proton/pion ratio 
at p t > 1 — 3GeV. This effect, the manifesta- 
tion of the collective radial flow, was suggested 
and searched for but not found in [3] using the 
min.bias ISR pp data long ago. It works spec- 
tacularly in AuAu collisions at RHIC and PbPb 
collisions at LHC. The same effect should be 
looked at in the (triggered) high multiplicity 
pA and pp events. 



B. The viscosity effect on the radial flow 

Here we continue to discuss the radial flow 
using Gubser's solution [T!5], now adding the 
viscosity effect. The equation for the reduced 
temperature T — e 1 / 4 using the combination of 
variables 



3 (l +5 2 )3 /2« +2g ^YT^f+g 2 H = (19) 
dg 

where the last term contains a new prameter 



V 

£ 3/4 



^ 4 / 1/4 

s 3 



(20) 



1 - t l + r 
2t 



(18) 



For rj/s — 0.134 we will use as representative 
number H — 0.33. 

The equation is solvable analytically in terns 
of certain hypergeometric functions (see below) 
or numerically. For AA collisions we find that 
the role of the viscous corrections is truly neg- 
ligible. The curves are the same within their 
plotted width. (This is, of course, well known 
from numerical studies in the literature during 
the last several years.) For the pA and pp cases 
one can see a difference between ideal and vis- 
cous , as we show in Figs. [4] through the re- 
duced temperature dependence T = T/t at cer- 
tain positions. The viscous effect is maximal 
at early times, and then the viscous and ideal 
curves meet. As expected, the viscous effects 
are very small at the fireball center r = 0, and 
become much more noticeable at its edge, see 
the r = 3 curve. In fact in this case the vis- 
cosity completely stops the cooling (decrease of 
the temperature) for a significant time, thus de- 
laying the freezeout. 

The main conclusion of this section is that 
a "realistic" viscosity of the sQGP is so small, 
that it provides a rather modest correction to 
the radial flow, even for the pA and pp collisions 
under consideration. 



C. High angular harmonics 

If the effects of order l/R are not negligible, 
they should be included. Keeping the first or- 
der gradient of the velocities leads to the cel- 
ebrated Navier-Stokes hydrodynamics. As one 
includes the second order corrections, one get 
other known approximations such as the Israel- 
Stewart approximation. Recently, using the 
AdS/CFT approach about a dozen of lowest or- 
der coefficients in the gradient expansion were 
identified with alternating signs. An approxi- 
mate PADE-like re-summation of these terms 
was suggested by Lublinsky and Shuryak [2Tj . 
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FIG. 5: (color online) Squared amplitude dissipa- 
tion factor Pj^ (as it appears in 2-particle correla- 
tors) for r//s = .134 as a function of the azimuthal 
harmonics m for AA (black) solid, pA (blue) dash 
and pp (red) dash-dot. 



perturbation with the wave vector k is 
8Ta V {t,k) ( 2-qk 2 t 



P k 



exp 



3s T 



(21) 



Since the scaling of the freeze out time is linear 
in R or tf « 2i?, and the wave vector k corre- 
sponds to the fireball circumference which is m 
times the wavelength, then 



FIG. 4: (color online) The temperature versus di- 
mensionless time t, for ideal hydrodynamics (solid) 
and viscous hydrodynamics with r//s — 0.132 
(dashed) lines. The upper pair of (red) curves are 
for pp, the lower (black) ones for pA collisions. The 
upper plot is for r = 1, the lower plot for r — 3. 



We will discuss the role of these higher order 
gradient corrections in section [nil 

The effects of viscosity are likely to damp 
more the higher angular flow moments, as first 
discussed by Staig and Shuryak [22] and recenty 
applied to a wide range of RHIC data [23 . The 
"viscous filter" for the amplitude of a sound 



2ttR = 



2tt 



inserting these values in (21) yields 



Pn 



exp 



G) 



l 

TR 



(22) 



(23) 



The exponent contains the product of two small 
factors, -q/s and 1/TR times the harmonics 
number squared. For PbPb LHC collisions one 
finds that 



1 

TR 



0(1/10) 



(24) 



and a comparable r\j s. For the most central 
bin - so that the elliptic flow does not obscure 
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genuine higher harmonics via nonlinear effects 
- one can immediately see from this expression 
why we observed harmonics for to = 1..6, with 
m > 7 still in the statistical noise (see e.g. the 
comparison to hydrodynamics in |24| ) . 

Proceeding to pp and pA collision (at LHC) 
by keeping a similar initial temperature Tj ~ 
400 MeV ~ 1/(0.5 /to) but a smaller size R, 
results in a macro-to-micro parameter that is 
no longer small, or 1/TR ~ 0.5, 1, respectively. 
For a usual liquid/gas, with rj/s > 1, there 
would not be any small parameter left and one 
would have to conclude that hydrodynamics is 
inapplicable for such a small system. How- 
ever, since the quark-gluon plasma is an ex- 
ceptionally good liquid with a very small rj/s, 
which would keep this damping under control 
for to = 0,1,2, while harmonics with m > 3 
would be hard to see. 



nonzero viscosity the solution is 

f H sinh 3 p 



T : 



(cosh/?) 2 / 3 9(coshp) 2 / 3 

2 



X - F] 1 2^2 ;_Si ; ,27i 



with f = rfl /4 T and /* = e/T 4 = 11 as in 
PI- 

Small perturbations to Gubsers flow obey lin- 
earized equations which have also been derived 
in (2D]. We start with the zero viscosity case, 
so that the background temperature (now to be 
called Tq) will be given by just the first term in 
(27). The perturbations over the previous so- 



lution arc defined by 



T 



Ml + 6) 



(28) 
(29) 



D. Angular harmonics of Gubser flow 

Unfortunately, the formulae we have used so 
far are, strictly speaking, only true for asymp- 
totically high harmonics, with to 3> 1. Since we 
are actually interested in not so large to — 2, 3, 
we may check how the viscous filter works using 
a better approximation. For that, we return to 
Gubser's flow and consider its angular pertur- 
bations. The latters have been considered in 
the literature [20 1 124 ] . 

In [20] Gubser and Yarom re-derived the ra- 
dial solution by going into the co-moving frame 
via a coordinate transformation from the r, r to 
a new set p, given by: 



sinh p 
tan 9 



1 



2r 



2r 



1 



(25) 
(26) 



In the new coordinates the rescaled metric 
reads: 



ds 2 = -dp 2 + cosh 2 p ( 



sin 



6d(j) 2 ) + drf 



and we will use p as the "new time" coordinate 
and 9 as a new "space" coordinate. In the new 
coordinates the fluid is at rest, so the velocity 
field has only nonzero u p . The temperature is 
now dependent only on the new time p. For 



with 

u „ = (-1,0,0,0) (30) 

u lt2 = (0,ug(p, 9, (p), u ( j,(p, 6, (f>), 0) (31) 

5 = 5(p,e,4>) (32) 

The exact solution can be found by us- 
ing the separation of variables S(p, 9, <j>) = 
R(p)®(9)$(9). In the non-viscous case, that 
we are now discussing, each of the three equa- 
tions 

R(p) + \ tanh P R{p) + - 2 R(p) = 
a 3 cosh p 



M0) 



" ^ ^ tan u 

$(0) +TO 2 $(0) = 



A 



sin 2 9 



0(9) = 
(33) 



are analytically solvable, with the results dis- 
cussed in 24 . The parts of the solution de- 
pending on 9 and 4> can be combined in order to 
form spherical harmonics Yi m (9,(f>), such that 
S{p,9,4>)ocR l (p)Y lm (e, < t>). 

The basic equations for the p-dependent part 
of the perturbation, now with viscosity terms, 
can be written as a system of coupled first- 
order equations [20] • We are assuming rapid- 
ity independence, thus the system of equations 
(107), (108) and (109), from the referred pa- 
per, becomes two coupled equations, for (the 
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p-dependent part of) the temperature and ve- 
locity perturbations 

where the index v stands for viscous and the 
matrix components are, 



H tanh 2 p 
1(1 + 1) 



3Tt, cosh p 
2H Q tanh p 



, (^H tanh p — f b 



Tn = 
Ti 2 = 

r 2 i = 

Hq tanh p — 2T b 
T 22 = (8T 2 tanhp 

V cosh p 
+6iTo tanh 3 p)/(6T b (h tanh p - 2T b ) ) 

(35) 

Before we display the solutions, we need to 
translate our space-time plot into the p — 9 co- 
ordinates. The initiation surface t = ij are not 
the p = const surfaces. The freezeout ones also 
do not correspond to fixed p because the tem- 
perature is T = T(p)/t(p,9). So, in both cases 
one has to decide which points on the initiation 
and final surfaces are most important. The thin 
solid lines in Fig[3] approximately represent the 
initial pi and the final p/ values for all three 
systems. Therefore, we will solve the equations 
between those two surfaces. 

In Fig|6] we show the solution of the p evo- 
lution of the two variables, the temperature 
perturbation and velocity 5i (p) , vi (p) . As one 
can see, all of them start at p = — 2 from 
the same Si = 1 value. While the elliptic 
one I = 2 (black solid curves) changes more 
slowly, higher harmonics oscillate more. The 
ratios v m /e m are predicted by the hydrodynam- 
ical solution, complemented by the Cooper-Fry 
freezeout. Since the latter is rather involved 
(see [24] for discussion in detail) we will rea- 
son on the basis of the former alone. Assuming 
that the perturbation amplitudes are relatively 
small and everything scales approximately lin- 
ear, and vi at freeze out dominate, one can re- 
late the relative magnitude of the harmonics 



into the measured two-particle correlation func- 
tions to scale as squares of the flow harmonics 



(v 



K A ) 2 



(v. 



om4 a ) 2 



0.3(ef; 



PP-yi 



2 . 



: 0.16(ef Y (36) 



This is qualitatively consistent with the 
squared damping of the amplitude of the pre- 
vious section, for m = 2. 

A comparison to CMS data shows that the pp 
data show smaller vi as compared to pA data. 
Quantitatively, the ratio is about a factor of 1/4 
(see Fig. 3 of [7]) rather than 1/2 which the hy- 
dro solution provides. Perhaps it is because the 
pp collisions create a somewhat more spherical 
fireball, with e 2 p < , in spite of having a 
smaller size. We will return to this issue at the 
end of the paper. 

Let us now compare in a similar manner the 
ratio of the m = 3 to m = 2 harmonics 



,AA 



,,AA 



0.12 



0.09 



0.02 




(37) 



Assuming €3/62 ~ 1 one finds that in pA we 
predict W3/W2 ~ 1/3, which agrees nicely with 
the ALICE data |S]. For pp we have v^/v 2 ~ 
1/7 which is probably too small to be seen. 

One may further ask why the experiment 
shows a very large 1*3 jvi > 1 in AA? The reason 
is explained in detail in |24) . It is related with 
the issue of the "sound horizon, and a larger 
temperature correction in the third harmonics 
in comparison to the elliptic term 



Sj A (Pf) 
5£ A (po) 

S AA (Po)J 



0.1 



(38) 



which results into a very large contribution at 
freeze-out. This feature pointed out in [21] is 
nicely confirmed by the AA data (but not all 
hydro works) . As one can see from our plots of 
5i(p), it is not the case for the pA and pp cases. 
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AA 




P P 

FIG. 6: (color online) The dimensional less temperature perturbation 5i (p) and velocity vi (p) , for / = 2,3,4 
shown by (black) solid, (blue) dashed and (red) dotted curves, respectively. Three sets of calculations 
corresponds to A A, pA and pp collisions. 



III. HIGHER GRADIENTS 

A. The LS resummation 

The Navier-Stokes approximation used so far 
only includes the first order terms in the gradi- 
ent expansion. In this section we qualitatively 
discuss the role of the higher derivatives. Let 
us define symbolically the corresponding con- 
tribution to the stress tensor as 



n n ^ ' 

with coefficients c n and certain kinematical 
structures with i derivatives P£ v ■ Their order 
of magnitude is given by the pertinent powers of 
the small hydro parameter 1/TR. When these 
contributions are not small, they should be as- 
sessed, since they may cause the expansion to 
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diverge. 

AdS/CFT is an indispensable tool - in fact 
the only one we have - to address this issue 
for the sQGP. A consistent procedure deriv- 
ing hydrodynamics and including in principle 
any number of gradients from Einsteins equa- 
tion has been outlined in Ref. [26 . In prac- 
tice, it has been worked out to second order 
only. However for small (linearized) perturba- 
tions the correlators of the two stress tensors 
can be calculated to higher order in frequency 
and wave vector w, fc, extending the original cel- 
ebrated 77/s = 1/47T calculation [S3] to about a 
dozen further coefficients. 

An approximate PADE-likc rc-summation of 
the higher order terms has been suggested by 
Lublinsky and Shuryak (LS) |21j . The main 
point is to notice the alternating signs of the 
series, which calls for an re-summation a the 
geometrical series. Here we discuss only the 
single pole resummation model or LS2 in |21) 
in which the Navier-Stokes viscosity or NS is 
subtituted by an effective viscosity 



Vns 

1 - r /2 , fc 2 /(27rT) 2 - iuj mA /(2TrT) 

(40) 



Note that ( 40 ) involves only two dimensionless 
coefficients 

m,o = -\ %,! = 2 - ln2 = 1.30 (41) 

approximately reproducing all known terms 
as well as the large-fc, lu behavior. The 
re-summation into the denominator suggests 
reduced viscous effect as k grows. As a result, 
the LS prediction is that in pp collisions one 
gets effectively a smaller viscosity than in AA! 

This conclusion may sound too good to 
be true, and recently one of us has studied 
the "strong shock wave" problem [37] in the 
AdS/CFT setting, using first principles (solv- 
ing Einstein equations) and comparing to the 
LS resummation scheme. While this problem 
is also "hydro-at-its-edge" type, the gradients 
of a shock profile have no small parameter 
k/2nT ~ O(l). The deviations between the NS 
and the exact (variational) solution of the corre- 
sponding Einstein equations were indeed found 
to be on the level of few percents only. Studies 
of time-dependent collisions in bulk AdS/CFT 



have found that the first-principle solution ap- 
proaches the NS solution early on and quite ac- 
curately, at the time when the higher gradients 
by themselves are not small [28 . 
Changing k 2 ,ui into derivatives 



iui/q —} 



( — ) "I 

dr r dr 

d_ 
dt 



(42) 



makes the re-summed factor (with the denom- 
inator) an integral operator. For any function 
of the coordinates f(t,r) we define the "LS op- 
erator" acting on a function / as 



OZs(/)= 1 



fd 2 f , 19/\1 



2(2ttT) 2 \dr 2 r dr J f 

+ (2-ln2)-^-?/-- (43) 
V ' 2ttT dt f V ' 

Schematically the resummed hydro equations 
look as 

(Euler) = rjO L s(Navier - Stokes) (44) 

where 0^5 is an integral operator. However, 
one can act with its inverse on the hydrody- 
namical equation as a whole, acting on the Eu- 
ler part but canceling it in the viscous term 

Oj)^{Euler) = rj(Navier — Stokes) (45) 

These are the equations of the LS hydrodynam- 
ics. Obviously they have two extra derivatives 
and thus need more initial conditions for solu- 
tion. 

Instead of solving these equations, we will 
simply check the magnitude of the corrections 
appearing in the l.h.s due to the action by the 
LS differential operator on the (ideal Gubser) 
solution used as a zeroth-order starting point. 
As one can see, large systems have a small 
q/T ~ 1/RT parameter and so these correc- 
tions are parametrically small. The issue is 
what happens "on the hydro edge", when the 
corrections have no formal small parameter. 

In Fig{7] we show the (inverse) action of (431 
on the zeroth other temperature profile of the 
Gubser flow as a function of r. We have used 
the freeze-out temperature Tf — 150 MeV and 
the indicated respective freeze-out times for pp, 
pA and AA. The higher gradient corrections for 
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forces any extra derivative to carry a particle 
momentum. As a result, the expansion param- 
eter of the n-th term is 



FIG. 7: (color online) The action of the LS operator 
Ols (43 1 on the zeroth order (non- viscous) temper- 

The three lines 



ature profile, the first term of \27\ 
correspond to AA (black) solid, pA (blue) dashed 
and pp (red) dash-dot. 



AA and pA are inside the few percent range 
from 1, while in the pp case the correction is 
larger, yet still in the 15 percent range. We thus 
conclude, that if the LS resummation repre- 
sents the role of the higher gradients, the overall 
corrections remain manageable even for the pp 
case. 



B. The momentum dependence of the 
harmonics and viscosity at freeze out 

The issue of higher gradients should also be 
addressed at the kinetic level, as emphasized 
by Teaney |29j . Equilibrium local distribution 
function of (quasi) particles should be corrected 
by the gradient expansion terms 



f(p)= /o(p) + /i(p)VK3,«„ 
+ (higher gradients) 



(46) 



which are negligible in the macroscopic limit 
TR — > oo but not in the "hydro-at-its-edge" 
problems like the pp case we discuss here. The 
main point in [29 is that Lorentz covariance 



1 ( p 1 
s\TTR 



(47) 



In AA collisions the smallness of the second fac- 
tor allows p/T to become large or O(10), before 
the macroscopic theory breaks down. Indeed, 
the higher harmonics of the flow agree with this 
estimate for transverse momenta of the order of 
Pt «3GeV, or p t /T f « 20. 

Following this line of reasoning, and turning 
to the case of pp and pA collisions, the sec- 
ond parameter is no longer small. Therefore we 
expect the validity region of the macroscopic 
theory to be strongly reduced, say to a much 
smaller region in pt/Tf ~ 0(1). 

We know it cannot be correct, as both CMS 
and ATLAS are not even able to observe soft 
particles. In fact, as seen from Fig.4e of [5], the 
elliptic flow V2{pt) in pA rises linearly to about 
2 GeV/c, where the presummed viscosity effects 
cause it to start decreasing. Furthermore, the 
elliptic flow reaches there a magnitude compa- 
rable with that in AA collisions. So, apparently 
the experiment does not see an expected reduc- 
tion of the hydro validity region! 

The only way we think the data can be con- 
sistent with theory is because the viscosity rj/s 
should not be treated as a universal constant. 
While its value extracted from hydro calcula- 
tions of v n is an average over the whole duration 
of the collision, the viscosity we discuss now is 
at freeze-out. Perhaps its effective value there 
is even smaller than the average. 



IV. FROM THE USUAL TO 
EXPLOSIVE POMERON 

A. The SZ model: classical strings and 
Polyakov tachyon contributions 

The SZ model [HE] is based on bosonic string 
exchanges between the colliding high energy 
objects. It is essential that the QCD string 
with a nonzero tension related to QCD con- 
finement is used, and not the conformal super- 
string which has a massless spin-2 graviton ex- 
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A 



(a) 



b 

(b) 



FIG. 8: Dipole-dipole scattering with separation b: 
pomeron exchange (a); reggeon exchange (b). 



citation. There is no supersymmetry and gravi- 
tons transmutes to a massive spin-2 glueball 
with an exponentially small contribution in the 
pomeron diffusive limit [14j [17] . However there 
is still a large N c parameter, related with a 
small string coupling g s and a large 't Hooft 
coupling A = g s N c so that 1/A effects of the 
curved geometry will be considered as sublead- 
ing. 

At very high energies the rapidity interval pa- 
rameter can be used as a large parameter 



X = ln(s/s ) > 1 



(48) 



It will play the role of the effective time in what 
follows. Transverse momentum transfer is held 
fixed t = —q 2 and soft. The main phenomenon 
to be studied is the string diffusion. Two lon- 
gitudinal directions - time and the beam di- 
rection, also often used as light cone variables 
x± - are complemented by two transverse co- 
ordinates plus a "scale coordinate" z. Its ini- 
tial value corresponds to a physical size of the 
colliding dipoles and diffusion means the pro- 
duction of small size closed strings. The z- 
coordinate is not flat. We will model its metric 
by an AdSs with a wall. The number of trans- 
verse coordinates, which will play an important 
role in the following, is thus 



D, = 3 



(49) 



We will now review the pomeron results and 
its associated entropy in this setting. The am- 
plitude of the elastic dipole-dipole scattering in 
Fig. [§i reads [SHQ1] 



-2is 



T(s, t; k) « g 2 s J d 2 b e tq b K T (/3, b; fc)(50) 



where Kj- is the pomeron propagator for dipole 
sources of color A^-ality k describing the string 
flux, k runs over all integers till N c /2 for even 
N c and N c /2 + 1/2 for odd ones. In the real 
world with the SU(3) color group, k = 1 is 
the usual string between fundamental charges 
(quarks) and the largest tension k = 2 is the one 
between two baryon junctions. The first argu- 
ment of the propagator is f3 — 2ith/x, where 
b is the impact parameter. g s w 1/N C is the 
string coupling. 

The explicit form of Kj- for the standard long 
strings regime 

b > P > (3 H (51) 
follows from the Polyakov string action, 



K 



,b;fc) 



P 



D ± /2 



47T 2 b 

d( n ) e -' T ' 9b ( 1 -' 3 H/2/9 2 +47rn/^/3 2 ) 



(52) 



n=0 



We have defined a = or/2, the Hagedorn in- 
verse temperature as 



and the inverse temperature as 

27rb 



Pv = P/k 



xk 



(53) 



(54) 



The subscript U is a reminder that it is basi- 
cally the Unruh temperature related to a fixed 
acceleration/tension. The asymptotic string 
density of states is 



d(n) w e 2 -V^W^ /n D ± /4 



(55) 



with the normalization d(0) = 1. 

For high energies and large b - long strings 
- the pomeron propagator Kj- to leading order 
in 1/A obeys a diffusion equation in rapidity x 
and curved transverse space 



(d x + *D k (Mg-V£))K T = 



(56) 
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with Dfc = a'/k = l 2 s /k the pomeron diffusion 
constant. This diffusion (561 is reminiscent of 
Gribov diffusion and implies on average (b 2 ) = 
DfcX for close pomeron strings. For large b and 
confining AdS backgrounds, the diffusion takes 
place near flat geometry (BH horizon or wall) 
with a tachyon mass related to string modes 



4D ± 



' OO 



\n=l 



i 



(57) 



with D± = 3 in AdS$ with a wall. While we 
describe quantum string dynamics at zero tem- 
perature, the formulae resemble thermodynam- 
ical ones. This happens because the string- 
membrane configuration in question, depicted 
in Fig|8ja) is a tube, and quantization on a cir- 
cle is formally identical to thermal Matsubara- 
time formalism. 

The extra z coordinate is different from oth- 
ers. A finite size dipole sitting at a height z a 
finite distance from the confining wall, experi- 
ences corrections 0(A -1 / 2 ) due to the curvature 
in z. This results in a shift of the tachyon mass 



Mi 



(D± I) 2 
4a' VX 



(58) 



Most of the arguments to follow will be carried 
out for large A> 1 unless indicated otherwise, 
so this effect is considered small. 

Inserting the leading n = contribution of 



( 52 1 in ( 50 1 yields the pomeron contribution to 



the elastic dipole-dipole scattering amplitude at 
large x an d fixed N-ality k 



T(s, t; k) rj ig 2 s 



•so 



(59) 



Thus the resulting pomeron has the intercept 
above 1 (and corresponds to cross section grow- 
ing with energy) 



ap, fc (0) = 1 + 



kD ± 



1 



fcL> ± 
12 



(D± I) 2 
12 ' * ' 12 sV\ 

(60) 

where the first term is due to Casimir-Luscher 
contribution and the 1/y/X correction follows 
from the tachyonic correction |58| as discussed 
in®. 



The subcritical string regime discussed so far 
is defined by the condition f3 = 2nh/x > Ph 
which translates to a rapidity range (collision 
energies) x < 2 in the diffusive limit (fe 2 ) = 
-DfcX- A more precise bound follows from the 
inclusion of the 1/A corrections in the tachyon 
mass (58) or 



p > ^2(«p ~ 1) Ph 



(61) 



This leads to the bound x < 10 for the cor- 
rected phenomenological value of the pomeron 
intercept ap — 1 = 0.08 in (60), which roughly 
corresponds to energies below the LHC. This 
condition discriminates between a sub-critical 
and a critical string as we will detail below. 
We note that (61) implies a strong coupling 
rcnormalization of the Hagedorn temperature 
through the geometry of AdS. 



B. Reggeon 

For completeness, we note that reggeon ex- 
change with open strings can be addressed sim- 
ilarly. For the Reggeon a = a? and the elastic 
scattering amplitude for dipoles of N-ality k as 
shown in Fig. [8] b is now 



T(s,t;k)Kig 2 s f 50 ) (- 
\ s J V s 



(62) 



with the extra so/s pre-factor accounting for 
the normalization of the spinors traveling on 
the exchanged world-sheet. This point was 
originally made in [TJ] but with different con- 
clusions for the reggeon intercept. At large 
s, the pomeron exchange is dominant. The 
pomeron as a closed string can be viewed as 
2 glued open strings or a pair of reggeons up 
to spin factors. As a result the reggeon slope is 
twice the pomeron slope while its intercept is 
also twice the pomeron intercept. 



C. Entropy 

The pomeron described here is a non-critical 
string in D± +2 = 5 dimensions in leading or- 
der in 1/A. We already commented how the 
temperature is formally introduced: the same 
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can be done for the free energy F = —IioKt/Pu 
and the entropy [TJ [5] 



n=l ^ 

'(3 k 



+D ± 



12 2 V V 2tt 



1, 
(63) 



At large rapidity x the entropy is dominated by 
the tachyon contribution 



D±(3 k 
12 



(64) 



Since (3 k — the entropy scales with the 

rapidity interval x- I n contrast, the energy E ps 
ah with on average (b 2 ) ps DfeX: scales with 
the root of the rapidity interval, and therefore 
is subleading for asymptotically large x- This 
is a major difference between this regime and 
the explosive string that we will discuss below. 

Large particle multiplicities in pp, pA and 
AA collisions based on (64) can be achieved 



through multiple pomeron exchanges [TJ |5] as 
illustrated in Fig. [9j The multiple exchanges 
start to interact as they diffuse transversely. 
Their number density in transverse space is set 
by the squared stringy saturation scale [JJ |2] ■ 
Although each exchange is penalized by g 2 s ps 
l/N%, it yields to a shadowing of the dipole- 
dipole cross section and saturation after an 
cikonalized resummation. The multiplicities 
are found to be consistent with the bulk multi- 
plicities at current collider energies. 

In this scenario, the exchanged strings are 
sub-critical with T < Th- They may turn into 
an ideal fluid of string bits with very small 
viscosity [2]. The stringy matter is released 
sub-critical and non-explosive. It can trans- 
mute to a black-hole on a time scale of order 
t ps (b/x) 3 /(4Z 2 ) after the percolation in trans- 
verse space p. 



D. Near-Hagedorn strings and the 
Nambu-Goto description 

We recall that all the expressions discussed 
above such as (52) for instance, were derived 




FIG. 9: Multi-pomeron contribution to dipole- 
dipole scattering near saturation with sub-critical 
string exchanges. 



regime (3h < (3 < b dominated by the tachion 
mode n = 



K T (Ab;fe) 



47T 2 b 



D ± /2 

(65) 

However, as the Unruh temperature becomes 
closer to the Hagedorn temperature or (3u — > 
(3h, string excitation is no longer small. Thus 
the use of the Polyakov action is no longer valid 
and (3h/(3 corrections need to be resummed. 
For D± = 3 we have (3h — 2ttI s so that 
(3h ?s Pu requires b ps xh/k. This shortcoming 
is readily fixed by noting that ( [65] ) is just the 
large (3 asymptotic of the Nambu-Goto (NG) 
string action. The NG string is a conformal re- 
alization of the Polyakov string which is valid 
for smaller b and (3. The re-summed result ob- 
tains a square root (which will play important 
role in what follows) 



K T (Ab;fe) 



47T 2 b 



D ± /2 



2 /o2\l/2 



-a/3h(l-/3 2 H /p 2 ) 



using the scalar string or Polyakov action in the 



(66) 

In the NG realization of the pomeron the tachy- 
onic contribution to the transverse propagator 



17 



is still dominated by the tachyon provided 
that the impact parameter is larger than the 
critical be = ttI s . Clearly ( 66 1 reduces to (651 
for fa/0 < 1. 



The resummed expression ( 66 ) shows that 



the Unruh temperature causes the string ten- 
sion to vanish at the Hagedorn point 



a(l-/3|//? 2 ) 1/2 ^0 



(67) 



in agreement with the universal behavior ob- 
served for static thermal strings [3T]. As we 
noted above, this occurs when the impact pa- 
rameter b w xls- 

The scattering amplitude associated to the 
NG tachyon can be obtained by inserting ( 66 1 
in ( 50 1 . The result in the saddle point approx- 



imation reads 



T(s,t; 1) w ig s 

(t/v / 2)(l-i(l + ^/l-2/t)) (l+^/l-2/t+l/t) 



SO 




RH 



FIG. 10: The string ball: RH is the Rindler horizon 
along the longitudinal direction L, with the string 
streching along the transverse T direction with im- 
pact parameter b. 



with t = a't and k — 1. (69) reduces to the 
pomeron amplitude (59) for s 3> — t > 1/a'. 



The cross section ohm for the production of 
near-critical strings or high multiplicity events 
follows from ( 50 1 , by the optical theorem , with 
q = and b w %i s . Specifically, for dipole- 
dipole scattering with N-ality k = 1, one hnds 



ohm 



( Ph 



Cl/2 

Ph 



(69) 



with 



A/3 



>H 



Th 
T 



1 = O 



(70) 



The minimum bias cross section gmb has been 
estimated in [17] after an eikonal resummation 
of the subcritical strings, 



TtD^a' 2 

&MB ~ — n — X 



(71) 



The ratio of the high multiplicity events to the 
minimum bias events can be estimated as 



~9s 



D ± /2 



X 



A/3 
Ph 



(72) 



where we have dropped an overall number of 
order 1 and D± = 3. For near critical strings 
we have ([70} and g s w l/N c . At LHC, x ~ 10 
so that 



ohm 
a mb 



10" 



(73) 



This estimate is comparable to the probability 
of the high multiplicity trigger used by the CMS 
collaboration for events displaying the ridge in 
pp collisions at LHC. 



E. The string ball as a black hole 

We now show that the string at this point 
transmutes to a string ball (dual to a black- 
hole) with large energy and entropy. The near- 
Hagedorn pomeron has a horizon located along 
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the longitudinal L-direction shown in Fig. |10| 
This is very different from the effective ther- 
mal AdS metric used in most of the thermal 
AdS/CFT studies of the sQGP because the 
horizon is not along the holographic z-direction. 
We note that the SZ model has a holographic 
z-coordinate in the T-transverse coordinates in 
Fig. |10| but with a confining wall at z — zq, not 
a thermal horizon. 



The pomeron propagator ( 66 1 behaves like 



a thermal ensemble with Unruh temperature 
1/Pu- We follow pQ and define the free energy 
or pressure F = — lnK-r/^u, 



F(/3,b) « kah 1- 



1/2 



(74) 



At the near-Hagedorn regime with b w xh/k 
the pressure is small, but the energy and en- 
tropy are very large 



E = dfaiPukF) « kah 



p 2 



(75) 



,2« v~r 6 l/0)ko-h(l-^ V2 



S = p 2 d 0u F « 



For /3 rj the first law of thermodynamics 
for black-holes in Rindler coordinates as noted 
by Susskind [32|, emerges with 



(76) 



and vanishingly small pressure (74 1. We note 



the Rindler temperature Tr = 1 /2tt and there- 
fore the Rindler energy E# = EZ S . The emer- 
gence of a Rindler temperature is expected 
since the stringy pomeron exchange is charac- 
terized by a line element [T] 



ds 2 



a 2 p 2 di A + dp z 



dsi 



(77) 



with a Rindler acceleration a = x/b. At this 
regime the acceleration is a — k/l s . On the 
streched horizon at p = l s /k in (77), the warp- 
ing of time is 1 since t/t p = (b/x)/p —> l s /kp. 
A cartoon of the string ball as a black-hole is 
shown in Fig. [TO] 

The transverse area of the black-hole is the 
area of the diffusing string in rapidity 



.4 



BH 



(78) 



in transverse D± = 3 provided that the diffu- 
sion length in the z-direction is within the con- 
fining wall. As a result, we have the Bekenstein- 
Hawking type relation 



1 

AG~ 5 



A B h 

with an effective Newton constant 



G 5 



r 2 (( X A 3 )(l-/3f,//? 2 )) 1/2 ll 



(79) 



(80) 



For a fundamental string, the Planck and string 
constants tie with G 5 through G5 = lp = g 2 l^, 
with in general g 2 ps 1/N 2 . 

The transmutation of the fundamental string 
to a black-hole at the string scale was fore- 
seen by Susskind and others in the context of 
string-based gravity [321 133] ■ Indeed, through 
an isentropic variation of the string coupling 
g s — > as a trick, it was argued that the 
Bekenstein-Hawking formulae emerges from a 
statiscal counting of quantum string states. In 
hadronic collisions at large rapidity x, the effec- 
tive relation ( |80| shows that this transmutation 
can be achieved in a twofold way: 1/ (3/ Ph 1 
(Hagedorn limit); 2/ kj\ — > 00 (large N-ality 
or N c limit). The Hagedorn regime can be 
achieved through fluctuations with high mul- 
tiplicity events at current collider energies in 
pp, pA and AA experiments as we suggested 
above. 

Empirical estimates based on DIS data anal- 
ysis suggests zo ~ 2/GeV [T], so that the dif- 
fusion length is within the confining wall for 
\J x/kl s < zq or x < 16 for k = 1. For x > 16 
the diffusion length is outside the confining wall 
and the transverse area ( 78 ) is now changed to 



.4 



BH 



2tt 2 z 



(81) 



with the corresponding changes in the effective 
Newton constant estimate 



G, 



(n 2 /k)(l-p 2 H /p 2 ) 1/2 (z l 2 ) (82) 



F. String on the Rindler horizon 

At the near-Hagedorn regime the pomeron 
string carries large energy and entropy but van- 
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ishingly smal pressure. To assess the transverse 
viscosity of this string without recourse to a full 
Nambu-Goto analysis, we make the key obser- 
vation that in near Hagedorn regime the expo- 
nents in ( 63 ) transmute to 



f3 k n -> f3 R n 



(83) 



with n the Rindler frequency of the harmonic 
oscillator. At large l/fik or large Rindler tem- 
perature 1//3r, the string energy and entropy 
are dominated by the high modes in ( 63 1 . We 



have explicitly checked that the tachyon ther- 



modynamics ( 74 ) and ( 76 1 follow from the large 
n excitation spectrum of the NG string by us- 
ing the modular transformation and the saddle 
point approximation, as we detail in Appendix 
B. 

At the Hagedorn limit, a long and space fill- 
ing string, with D± dimensions, is a very effi- 
cient way to carry large entropy. The analogy 
between a string ball and black hole thermody- 
namics shows that in fact it carries the largest 
entropy density possible! With this in mind and 
for simplicity, consider a Polyakov string made 
of D± harmonic oscillators immersed in a heat 
bath with finite but large Rindler temperature 
1/ (3r. The energy of the string is dominated by 
the high frequency modes, 



Efl ps D ± 



71=1 



,J3n.n _ i 



(84) 



For large 1//3r it is black-body 



E, 



(85) 



Through the first law of thermodynamics ( 76 1 



we can enforce the zero pressure condition on 
this highly excited string, with 



TT 2 D A 

2^T 



(86) 



G. Viscosity at the Rindler horizon 

Viscosity can be defined via certain limits of 
the correlators of the stress tensor, known as 



the Kubo formula. Thus one does not need 
hydrodynamics to calculate it, just the stress 
tensor. To assess the primordial viscosity, we 
follow [2 and write the needed expression on 
the streched horizon for the excited string 



VR = lim 

UR—yO Zojr 



dre 1 ' 



R 
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(t) (87) 



with Ar the area of the black-hole and r a di- 
mensionless Rindler time. The retarded com- 
mutator of the normal ordered transverse stress 
tensor for the Polyakov string on the Rindler 
horizon reads 



R 23 ,2 3 (r) = ([Tf(r),Tf(0)]) 



with 



T?(t) 



2A, 



2 ^ 



-2im 



(89) 



and the canonical rules 



, a? n ] = mS m+nfl 6 ZJ . 
The averaging in (|88|) is carried using the black- 



body spectrum as in (84). The result is 



Ar 

rjR = lim - — 

IdR—tO ZUJR 



(W2) 5 



l 



1 



Ar 8(3r 
(90) 

We note the occurence of the Bekenstein- 
Hawking or Rindler temperature Pbh = Pr in 
the thermal factor. 



Combining ( 86 1 for the entropy to ( 90 1 yields 



the viscosity on the streched horizon 



Vn 



1 

4?r 



3 



l 

47T 



(91) 



which, for D± = 3, is precisely the celebrated 
universal value from AdS/CFT. The result (91 1 



is remarkable as it follows solely from a string 
moving at large "time" x m non-critical dimen- 
sions but near its Rindler horizon, not in trans- 
verse coordinate z. It emerges naturally in the 
near-Hagedorn regime. 

The result (91 1 for the critical pomcron as a 



close string exchange on the streched horizon 
for large 1/Pr is to be contrasted to the same 
viscosity ratio but for the low-T pomeron as a 
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close string exchange far from the horizon for 
small l//3 fc 



H. Very high multiplicity and transition 
to sQGP 



1 1 

24^ 



(92) 



The ratio is small at large rapidity. ( 92 1 re- 
duces to (91) for fik —> Pr up to a factor of 



1/2, showing the non-commutativity of the two 
limits. Indeed, for small l/fik the non-critical 
pomeron is described by the Polyakov action 
whereby the zero pressure condition (emblem- 
atic of a near-Hagedorn or black-hole in Rindler 
coordinates) does not hold. 



In light of (40) above it is amusing to note 



that the stringy relation (90) yields on the 



Rindler horizon an effective viscosity for finite 
frequency but zero momentum 



(93) 



( 93 1 thermally suppresses the large frequency 



modes ojr. The on-set of the black- hole is fol- 
lowed by black-body or Hawking radiation of 
string bits of frequency ujr/I as is explicit in 
(93) and stressed further below. In particu- 



lar, for finite wavenumber fc^ in Rindler units, 
the suppression is physically expected to follow 
from the substitution 



UR 



^R 



(94) 



and therefore exponential as well. The effec- 
tive viscosity ^(wij, /cr) at higher gradients - 
larger k^- would indeed imply a smaller effec- 
tive viscosity in pp than in AA, as we discussed 
earlier in the LS re-summation and also from 
phenomenology. 

Concluding this discussion of the viscosity let 
us point out that while we can use the Kubo for- 
mula for any setting in which the stress tensor 
is defined, the viscosity itself is of hardly any 
use outside the hydrodynamical setting. And 
alas, both for the cool subcritical strings and 
the near-Hagedorn critical strings we so far con- 
sidered, there is no hydrodynamics involved. 
However, this is going to change as we now dis- 
cuss. 



The typical (minimal bias) dipole-dipole - 
and thus pp, pA and A A collisions - involve the 
exchange of pomerons and reggeons in the sub- 
critical regime with low entropy /multiplicity 
and low T < Th- At very large ^/s the 
cross section and the string impact parameter b 
grow, and eventually the possibility of observ- 
ing string fluctuations in the near-Hagedorn 
regime becomes increasingly probable. 

Furthermore, as the entropy density grows 
further, the string eventually fills the space 
and its repulsive self-interaction becomes im- 
portant. In thermodynamics this is seen as the 
end of the Hagedorn regime (the high end of the 
mixed phase) from which the system goes into 
a deconfined phase. We now suggest that an 
individual string ball reaching such an entropy 
density (corresponding to the point B in Fig{2]) 
does the same. Its entropy density S/T 3 stops 
growing. From ( 79 ) the onset is expected when 



the (Bekenstein-Hawking) entropy becomes as 
large as that of the gluons, which coresponds to 



A 



BH 



1 

4Gs 



Nl 
0% 



(95) 



Note that at large N-ality k l/iV c , so that the 
transition takes place very close to the Hage- 
dorn temperature /3//3 H - 1 = 0(l/N c ). As 
soon as the string entropy gets so high, the pres- 
sure is no longer small, eventually reaching e/3 
as for a gluonic gas, so an explosion follows. 

We suggest that these explosive events are 
triggered in the current LHC collider in pp and 
pA events. The black-hole explosion is signaled 
by a large radial and elliptic flows in the large 
multiplicity events. 

Let us now estimate the critical multiplic- 
ity above which such explosion should happen. 
The multiplicity N associated with the explo- 
sion of the black-hole is readily estimated from 
the entropy near the Hagedorn transition sub- 
ject to the condition (95). Indeed, 



N « 7.5 S « cr p H h (1 - 0% IP 2 



-1/2 



(96) 



where the conversion factor of 7.5 is borrowed 
from the entropy-to-hadron density relation at 
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freeze out. Substituting Ph/ @ — 1 ~ 1 /A^ c from 
(95), we find the critical charge particle multi- 



plicity to be 



N«7.5S«7.5x/2«45 



(97) 



after using that b ss x^s and x ~ 12 at the 
LHC. The estimate ( |97| ) can be compared with 
the measured threshold multiplicity of N > 50 
for events with the ridge, according to the CMS 
collaboration [Sj . However, we note that the ac- 
tual acceptance of the CMS detector in rapidity 
and more importantly in p t , will increase the 
latter number by about a factor of 2 or so. The 
freeze out conditions are probably also different 
in pp and pA from AA (from which the 7.5 fac- 
tor is taken). Only further studies of spectra 
in "explosive" events will make a comparison 
quantitative. 



I. Hawking radiation 

The typical pp,pA collisions create strings 
which decay via Schwinger pair-production 
mechanism, which is described by the Lund 
model, PITHIA etc. The proposed black-holc- 
like description -dual to a string ball - we as- 
cribed to high-multiplicity events leads natu- 
rally to a different particle emission mechanism, 
namely via the Hawking radiation. This emis- 
sion is fully thermal from the prompt emission, 
and not delayed as in string breaking mecha- 
nism. Let us note that our approach is similar 
in spirit but different in details to the Unruh- 
Hawking effect discussed in [34 . In particular, 
in this paper the Unruh temperature is related 
to acceleration created by a string, not a col- 
lective "effective gravity" of the string ball as a 
whole, as in our model. 

The power spectrum or Hawking emission 
per unit time from a black-hole is generic. For 
our rapidly moving string it involves a black- 
hole in 1 + 4 dimensions with the extra dimen- 
sion accounting for changes in the dipole scales. 
In D± + 2 dimensions it reads 



d^ +1 P = 



(98) 



d D ^ +1 k 



eU /T B n + (-l)2s+l (2,r) 



D ± + l 



We have only kept the dominant S-wave con- 
tributions. Here T BH = T R = l/(2nl s ). The 
sum runs over the spin s of the emitted particle 
with a s (uj) the S-wave absorption cross section 
or grey-body factor of a spin-s on a black-hole. 
For ul s < 1 



BH 



i Ks lp ±+1 S 



(99) 



with Abh the area of the black-hole. The last 
identity follows from the Bekenstein-Hawking 
type relation and shows that the power spec- 
trum is extensive with the entropy. For scalars 
Ks = 1 |35) . As the Hawking emission through 



(99) unfolds, the mass and radius of the black- 



hole decreases, causing the Hawking tempera- 
ture Tbh to increase. The emission process is 
inherently a non-equilibrium one. Here and for 
simplicity we assume it to be quasi-adiabatic 
with (99) adjusting to the change in T BH . 
For massless particles u — \k\ in (99). The 



luminosity defined as L(o>) = dP/duj for D± = 
3 is 



L(w) 
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BH 
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e ui/T BH _|_ ("_2)2s+ 



T (100) 



It is a black-body spectrum from a 5- 
dimensional space where the black-hole origi- 
nated from. 

As many of these black-holes are expected to 
be released in AA collisions they are the seeds 
of the primordial matter viewed as a collection 
of these tiny black holes. Primordial Hawking 
emission of partonic constituents as well as elec- 
tromagnetic radiation is what current heavy ion 
colliders are probing. We recall that for x < 16 
we have Abh ~ X 3 ^ 2 while for x > 16 we 
have Abh ~ X because of confinement in the 
holographic or conformal direction of the string. 
Therefore, we estimate the tresholds Nr(x) f° r 
the large multiplicity events with explosive hy- 
drodynamical flow to scale with beam rapidity 
as 



N T (xi) 

N T (X2) 



3/2 



(101) 



f° r Xi,2 < 16 and 1 for xi,2 > 16, irrespective 
of whether it is pp, pA and A A collisions. 
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V. SUMMARY AND DISCUSSION 

High multiplicity pp and pA collisions are 
very interesting systems which are supposed to 
display a transtion from a "micro" to "macro" 
dynamical regime. In this paper we tried to ex- 
plain how this transition works using the lan- 
guage of both macroscopic and microscopic the- 
ory. 

We started in the "macro-to-micro" limit, 
using the hydrodynamical framework. As we 
emphasized in the Introduction, the applicabil- 
ity of hydrodynamics to high energy collisions 
rests on two small parameters: (i) the micro-to- 
macro ratio 1/TR, (ii) the viscosity-to-entropy 
ratio 77/s. For central A A collisions, both are 
small or of order 0(1/10). For high multiplic- 
ity pA and pp collisions, the first parameter is 
no longer small 1/TR = 0(1), prompting us 
to ask which hydrodynamical predictions are 
preserved by the smallness of only the second 
parameter 77/s. 

After solving the hydrodynamical equations 
we found that the radial (axially symmetric) 
flow is little modified by viscosity and is in fact 
enhanced by higher transverse gradients. Thus 
our main prediction is an enhanced radial flow. 
Its signatures are well known [1] with specifi- 
cally a change in the observed p t spectra on the 
particle mass, or growing proton-to-pion-ratio 
with p t . The magnitude of the effect should be 
even larger than in A A collisions, see (17) for 



specific values of the maximal transverse veloc- 
ity calculated in the Gubser flow. 

However, higher harmonics are more penal- 
ized by viscous corrections. We used the generic 
high-harmonics expression from |22j and ob- 
tained explicit solutions for Gubser flow for 
m = 2, 3, 4 as shown in Figj6j We have found 
a small v^fv-i pa 1/3 ratio for pA in agreement 
with the reported ALICE data (to be compared 
to the value above 1 in central AA). The value 
of t>2 itself is also suppressed by viscosity, and 
the relative suppression we have found between 
the pp and pA collisions agree reasonably with 
the CMS data. 

Then we discussed the problem in the oppo- 
site "micro-to-macro" limit, in the framework 
of a holographic QCD string pomeron. We have 
reviewed the general arguments behind the SZ 
pomeron as an exchange of a non-critical string 



in curved AdSs with confinement. For long 
strings with /3 > the pomeron follows from 
the scalar Polyakov action for the slightly ex- 
cited string oscillators. The Luscher term gen- 
erates the intercept of the pomeron ( 59 1 , which 
for D± = 3 and, with a finite 1/A correction, 
yields its phenomenological value of 0.08. 

The fluctuations in "hotter" strings show 
that their effective Unruh temperature may be- 
come close to the Hagedorn temperature of the 
gluodynamics strings. Because of the Polyakov- 
Susskind phenomenon, string excitations are no 
longer small. The string shape is better de- 
scribed by a Nambu-Goto string instead of a 
scalar or Polyakov string. For such a large am- 
plitude of excitation a "string ball" forms, with 
a large energy and entropy, yet a small pres- 
sure. The ensuing dynamics closely resembles 
what takes place in equilibrium thermodynam- 
ics in the mixed phase of pure gauge theory. We 
have shown that it is dual to a miniature black- 
hole with a Rindler horizon along the longitu- 
dinal direction of the collision. This black-hole 
is 5 dimensional, with 3 transverse coordinates, 
2 spatial ones and 1 conformal z describing the 
scale evolution. The black-hole radius and area 
are set by the Gribov diffusion length, which 
grows with the collision energy x = m (s/so) as 
X 1 / 2 , and x 3 / 2 for \ < 16, respectively. For 
X > 16 the area growth is reduced to \ be- 
cause of confinement along the holographic or 
conformal direction of the string. 

Using the Kubo formula for string excita- 
tions, we found that on the stretched or Rindler 
horizon the shear viscosity to entropy ratio is 
precisely 1/47T, for D± =3. It is the same as for 
AdS/CFT black hole, in spite of the fact that 
these two black holes are very different. Ours is 
dynamical with a horizon normal to the longi- 
tudinal coordinate, while the AdS/CFT one is 
static, with a horizon normal to the transverse 
holographic z-direction. 

We have argued that when the Unruh tem- 
perature narrows even closer on the Hagedorn 
temperature, with T/T H - 1 = 0(1/N C ), the 
Bekenstein-Hawking relation for the stringy 
black-hole becomes partonic and the ball of 
strings turns into a sQGP fireball. As a re- 
sult, transition to deconfined phase unleashes 
a large pressure with p ss e/3 and the stringy 
black-hole explodes hydrodynamically, follow- 
ing the general scaling of viscous hydrodynam- 
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ics in small volumes, as we discussed above. 
Our estimates put the trigger for such high mul- 
tiplicity events at LHC at 10 -5 which is com- 
parable to the currently used trigger of 10~ 6 for 
pp collisions used by the LHC collaboration. 

Our estimates show that the treshold mul- 
tiplicity for the explosive pomerons is about 
N ~ 45 for pp collisions with a rapidity in- 
terval of x ~ 12- After outlying the suggested 
theoretical scenario, let us focus on its possible 
experimental tests. The most obvious of these 
tests consists in further investigating whether 
the "ridge" is indeed a collective phenomenon, 
rather than a correlation of few particles. It 
is by now well known how to carry this: one 
should compare the elliptic flow parameters ex- 
tracted from 2,4,6 particles. Collective phe- 
nomena induce approximately factorizable mul- 
tiparticle distributions, while dynamical corre- 
lations (like mini-jets) disappear in many-body 
correlators. Another obvious direction is to 
look for the "ridge" and related phenomena at 
RHIC energies. 

Our main proposal is to search for the radial 
flow, using spectra of identified particles, which 
ALICE, PHENIX and STAR can do in their pA 
sample, or perhaps even in their pp samples. If 
successful, one may use hydrodynamics back- 
ward in time and derive the fireball sizes R for 
the pp and pA systems at different multiplic- 
ity and collision energy, to be compared with 
our "ball of string" scenario. Perhaps one can 
even create string-based event generators and 
predict the elliptic fluctuations of such a ball 
and compare it to the t>2 observed as the (now 
famous) "double ridge" . 

Finally, we now discuss some important dif- 
ferences between our string based scenario and 
other models, and emphasize the differences 
in their microscopic dynamics. A rather com- 
prehensive description of the LHC pp and pA 
data has been done recently by Dusling and 
Venugopalan [36] using 1- and 2-pomeron ex- 
changes made of perturbative pomerons of the 
BFKL type. For typical collisions our string- 
based pomeron and their approach would per- 
haps agree, although our pomerons follow from 
strong coupling and large N c QCD. 

In contrast, the physics of small-probability 
and large-amplitude fluctuations we have dis- 
cussed in this work is fundamentally differ- 
ent. Weak coupling logic leads to saturation 



and "glasma", an ensemble of classical glu- 
onic fields with Gaussian fluctuations. The 
string approach we use is remarkably differ- 
ent: while the gluonic fields are also strong, 
they are confined into thin flux tubes, which 
have a very unusual dynamics and thermody- 
namics. At certain near-critical conditions (re- 
alized also for equilibrium gluodynamics as a 
"mixed phase") strings get "soft" and produce 
large energy/entropy clusters we called "string 
balls". Such clustering to very high density is 
in a way opposite to saturation. String clus- 
tering instead of saturation is natural in effec- 
tive gravity. Indeed, the thermodynamics of the 
string balls can indeed mimic that of the black 
holes. Note that these ideas, while foreseen in 
general terms in [34] . became only quantitative 
with the AdS/CFT correspondence and other 
studies of gauge-string-gravity dualities in the 
last decade. Its dynamical details still need a 
lot of work to be understood. One key feature 
of this approach is the completely different view 
on thermalization. Instead of rescattering of in- 
coming gluons, as done in weak coupling, it is 
instead a quasi-instantaneous [V Hawking radi- 
ation of new quarks and gluons induced by the 
effective "gravity" of the cluster near its hori- 
zon. 

Speaking about experimental tests of those 
divergent scenarios, let us first note the (some- 
what forgotten) issue of quark chemical equi- 
libration. We know from the structure func- 
tions that at high energies that the excited sys- 
tem produced should be mostly gluonic, with 
a small quark admixture. Perturbative logics 
lead to the so called "hot glue" scenario [37] in 
which quarks appear several fm/c later. Recent 
glasma models lead to similar conclusions. In 
contrast, the prompt black-hole formation and 
its ensuing Hawking type radiation should in- 
clude both quarks and gluons promptly. This 
point is supported indirectly by hydrodynam- 
ical applications: they work well, assuming 
prompt quark equilibration. 

The way to check this is through the pro- 
duction of photons and dileptons. The effective 
gravity and black holes are for strong interac- 
tions only, so photons and leptons are only pro- 
duced through a secondary radiation by black- 
body quarks. At RHIC the best place to tell 
the "hot glue" scenario from the full QGP one, 
is to look at intermediate mass dileptons, with 
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a mass range of 1 — 3 GcV. Unfortunately, the 
contamination of these leptons by charm de- 
cays remains an unresolved issue, which re- 
quires high-resolution vertex detectors. 

Another suggestion is to study the statistics 
of fluctuations. If pA data are to be interpreted 
as we suggest, it means that with the order of 
a few percent probability a "string ball" clus- 
ter is produced (and explodes into the double 
ridge). If so, in AA collisions the probability 
of this phenomenon should be 0(1), as we have 
hundreds of nucleons going through the near- 
center of another nucleus. While being only a 
small part of the fireball, the stringy cluster is a 
significant local perturbation of the initial con- 
ditions, resulting in sound waves and correla- 
tions of flow harmonics. This phenomenon ob- 
viously is not included in the standard Glauber 
theory of fluctuations, and perhaps can be re- 
vealed by flow harmonics correlations, induced 
by a localized source |24j . 
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Appendix A: The first harmonics and the 
recoil mass 



The first harmonics of the correlation func- 
tion has a special role since it is dominated not 
by flow per se but by the momentum conser- 
vation. It is usually considered to be just an 
uninteresting observable. For example no dis- 
cussion of this point is given in Refs. [7HS]- 

Borghini and collaborators [30] had calcu- 
lated its magnitude assuming that all the trans- 
verse momentum of the trigger is absorbed by 
all secondaries equally. This implies that the 
fireball as a whole moves with certain small ve- 
locity P/M where M is its total mass. The 
conditional probability in the 2-particle distri- 



bution is 



/(palPi) = /o(pa) 1 



2P1P2,X 



N(p 2 )(l+v 2 ) 



(Al) 

where the trigger , by definition, moves in the 
x direction, 



v 2 



(pI-pI) 

(Pi +Py) 



(A2) 



and the correction containing it appears here 
because in the denominator one has \Px) not 
the whole p t . 

First, let us argue that an assumption that 
all particles in the fireball can share the trig- 
ger momentum cannot possibly be true, since 
it contradicts hydrodynamical casuality. In- 
deed, in a hydrodynamical regime a perturba- 
tions cannot propagate faster than the speed of 
sound. Thus N in the formula above should 
be no larger than the number of particles in- 
side the sound horizon, from a trigger emission 
event. We then predict that in the magnitude 
of the vi one should use N recoi i/N to t < 1< 

One can estimate it as 



N 



recoil 



N t , 



Hi 2c s T f 
R 2 AY 



1/3 (A3) 



where the first factor is the transverse and the 
second longitudinal, AY" « 10 is the total ra- 
pidity interval between the beams. 

Let us now have a look at the experimental 
data. The fit for RHIC PHOBOS data gives 



V u 



-1.510" 



(A4) 



Substituting V2 ~ 0.2 and the mean square 
transvserse momenta, plus integrating it over 
momenta p 2 one gets the number of particles 
which are getting the "recoil" 



Pi < Pt > 
Via (Pf) (1 + v 2 ) 



500 



(A5) 



Indeed, that it is significantly smaller than 
the total multiplicity of the events in the cor- 
responding centrality class 10-20% of RHIC 
AuAu collisions N tot « 3000. 

Our second proposal is to look not at the av- 
erage but the distribution over the value P(vi) 
and interpret it in terms of the "recoil mass dis- 
tribution" . It is expected to be very different in 
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different models, and perhaps such data would 
have some discriminative power. 

For pA and pp collisions: systems much 
smaller than the central AA should have 
smaller recoil mass and larger recoil velocity, 
thus larger first harmonics. While extracting 
it from back-to-back jet component is not 
so easy, we urge experimentalists to measure 
the first harmonics as it contains important 
information. 



ah 2 



n s 



2n - 1 



> 1 



(B2) 



for which (Bl ) is to exponential accuracy 



K r 08,b;A;) 



(B3) 



Appendix B: Large n spectrum of string 
excitations 



in agreement with the tachyon result above. 
The string energy at the large n saddle point 



(B2) is 



It is interesting to note that the tachyon ther- 
modynamics can be alternatively derived by us- 
ing the modular transform of the transverse 
string propagator, (nearly) exchanging 6 f> |5 
and exchanging the close and open string de- 
scriptions. Indeed, the modular transform of 



( 52 1 can be cast as 



K T (/?,b;fc)« (Bl) 

OO 

jr d{n) e-^ h (i-bc 2 /b 2 +2W^)^ 2 

71=0 

1 /2 

with b c = ((7rDj_)(12cr)) = nl s and the den- 
sity of states (55 ). The NG form has been sub- 
sumed. ( |B1[ ) is seen to diverge for /} < /3#. 
The divergence is controlled by a large n sad- 
dle point, 



E « ah 



2irns 
ah 2 



ah 



(B4) 



and the corresponding entropy is 



lnd(ns) = 2tt 



D±n s _ £_l 
6 4 



lnn s (B5) 



which is seen to satisfy the zero pressure con- 
dition S « /3hE in leading order. They are the 
tachyonic energy and entropy in the Hagedorn 
limit discussed above. This is expected since 
the modular transform allows us to cross from 
the j3 < h regime of long and close strings, to 
the P > h of short and open strings. The two 
descriptions match at the border b w fi. 
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